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We report measurements of the two-magnon states in a dimerized an- 

tiferromagnetic chain material, copper nitrate (Cu(NOs)2 • 2.5D2O). Using 
inelastic neutron scattering, we have studied the one- and two-magnon ex- 
citation spectra in a large single crystal of this material. We compare this 
new data with perturbative expansions of the alternating Heisenberg chain 
and find good agreement with these calculations. The data may also show 
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evidence for the recently proposed 5 = 1 two-magnon bound state (Phys. 

Rev. B54, R9624 (1996)). 

PACS numbers: 75.10.Jm, 75.40.Gb, 78.70.Nx 
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Competition between hopping and binding effects in elementary excitations is a general 
feature of low dimensional hard core systems. Such effects are predicted to be manifest 
in the structure factor of multiparticle continua, accessible using neutron and Raman 
scattering, as bound modes and enhancement of continuum scattering |I[ [0 . Experimen- 
tal systems of particular interest in this respect are the S = 1/2 alternating Heisenberg 
chains (AHC) |l[ and the new spin ladder systems |§ which are predicted to have bound 
modes below the two-magnon excitation continua. To investigate this phenomenon we 
present an experimental study of the two-magnon states in a near ideal example of an 
AHC compound, using inelastic neutron scattering. 

The 5 = 1/2 AHC spin Hamiltonian is 

H — 22 J Sw-i ■ S 2 i + aJ S 2 i ■ SW+i (1) 

i 

where J > is the intradimer coupling, and aJ (0 < a < 1) the interdimer one, which 
alternate between chain sites i ||. Computational methods are very effective at calculating 
perturbative properties of (1), and ground and low-lying excited state wavefunctions are 
given to 0(a 5 ) in ||, together with many experimentally important quantities to 0(a 9 ). 
Higher order expansions for selected quantities have recently been reported Here we 
compare these results to Cu(NOs)2 • 2.5D 2 0, or CN for short. 

The magnetic properties of CN are well characterized. It is monoclinic (712/cl ||), 
with low temperature lattice parameters a = 16.1, b = 4.9, c = 15.8 A and (3 = 92.9° fl6[, 
with spin 5 = 1/2 magnetic Cu 2+ ions. The dominant magnetic exchange integral J is 
between pairs of spins forming dimers. Dimers separated by u are coupled together by 
exchanges J' . Only the J' a ,1 1, exchange paths are of appreciable strength, giving two 
sets of alternating Heisenberg chains (AHCs) running in the [111] and [111] directions 
which repeat every u = [lll]/2 and Uq = [lll]/2 respectively. 

Bulk magnetic measurements give information on the gap, exchange and ground state 
energies: Applied magnetic fields induce spin flop ordering in CN above a critical field 
H& 27 kOe, with a transition to full alignment at ~ 43 kOe. Because the orbital 
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moment is quenched, and demagnetization effects are negligible, H c \ directly gives the 
excitation gap, A = 0.378 ± .007 meV J7| and H c2 gives the sum of exchange couplings 
J + Eu = 0.580 ± .007 meV. Normally it is not possible to measure the ground state 
energy of quantum antiferromagnets but high field magnetization techniques in our case 
make this achievable. 

The ground state energy-per-spin eo can be found from the low temperature isothermal 
magnetization M(H) using eo ~ e/ — SH c2 + f" c2 M(H)dH where the fully aligned 
energy-per-spin is e f = S 2 /2 ■ ( J + £ u J' u ) = H c2 /8 = 0.0725 ± .001 meV. Using the 270 
mK data of Diederix et al. in Figure (3) of |7j (measured using proton resonance) to 
determine the integral over magnetization gives an experimental ground state energy-per- 
spin eo = —0.174 ± .004 meV. This is essentially the T = result, as the gap activation 
energy corresponds to 4.4 K. To estimate thermodynamic properties we approximate the 
sum of the exchanges by the single coupling aJ = Y^ U JL °f equation (1). Using the 
0{a 9 ) expansions § for A(a) and e gives J = 0.455 ± .002 meV and a = 0.277 ± .006; 
in agreement with the results of []7| and [|S|, J = 0.45 meV and a = 0.27. Our calculated 
values of the thermodynamic parameters J + J2u J' u = 0.581 meV, A = 0.379 meV and 
eo = —0.172 meV agree within error with the experimental values. 

The neutron scattering structure factor S(Q,cu) probes the ground \G) and excited 
states \E) of a magnetic system through the matrix element | (E\S a (Q)\G)\ 2 of the Fourier 
transformed spin operator S^Q) (a = x,y,z). We measured <S(Q, uj) using inelastic 
neutron scattering from two deuterated single crystals of CN with a total mass of 14.1 g, 
using the SPINS cold neutron triple-axis spectrometer at the NIST Center for Neutron 
Research. The substitution of D for H reduces incoherent scattering of neutrons and does 
not significantly change the magnetic properties of the material. The sample was mounted 
with (h, 0, /) as the scattering plane in a pumped 3 He cryostat at a base temperature of 
300 mK. 

The spectrometer was set up with 80' before the sample as the only collimation. A ver- 
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tically focused pyrolytic graphite PG(002) monochromator and a horizontally focused an- 
alyzer array composed of eleven independently rotatable PG(002) blades were employed. 
A cooled Be filter before the sample removed higher-order contamination from the beam. 
Measurements were made with fixed final energy Ef = 2.5 meV by scanning incident 
energy at various reduced wavevector transfers along the chain, q = Q ■ u . The wide 
angular acceptance (14°) of the analyzer dominated the instrumental resolution making it 
highly elongated perpendicular to the scattered wavevector in the scattering plane. Scan 
trajectories were chosen to maintain the final wavevector k/ along the (101) direction so 
as to integrate over nondispersive directions while maintaining good resolution in q. 

Because the ground state is a singlet it cannot be probed directly by neutrons, but its 
composition is reflected through spin matrix elements to triplet excited states. Neutron 
scattering matrix elements to the 5* = 1 one magnon states have been calculated to 0(a 5 ) 
0: the leading order scattering process is from the bare dimer component of the ground 
state, and an a/2 ■ cos(g) component in the one magnon structure factor arises from an 
0(a) two-dimer excitation in the ground state. Similarly, transitions between various 
components of the full ground and excited states can be identified through distinctive 
modulations of the Q dependence of <S(Q, uj). 

Figure (1) shows scans in energy in CN: Panel (a) shows a scan at the antiferromagnetic 
zone-center, q = 2tt, taken at T = 300 mK. Strong elastic scattering from incoherent 
nuclear processes is clearly seen as well as a one magnon || peak at 0.4 meV, close 
to the dimer energy J = 0.45 meV. The non-magnetic background (dashed line) was 
modelled by a Gaussian (incoherent) component and a power-times-Lorentzian (broad, 
quasielastic) component. A second magnetic peak appears at about 0.9 meV, double 
the dimer energy. Panels (b) and (c) show this peak with the nonmagnetic background 
subtracted at q = 2ir and q = 3ir, respectively. This feature is considerably weaker than 
the one-magnon scattering and narrows at the zone-boundary, panel (c); this behavior is 
consistent with two-magnon scattering. 
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Figure (2) shows excitations calculated to 0(a) using degenerate perturbation theory, 
including states up to two-magnon, for J = 0.45 meV and a = 0.27. The calculation 
is similar to that for the two-soliton continuum in the XXZ Ising chain in fl9|, except 



that a dimer excited state basis is used and the one magnon band is included [I0|. An 
interesting feature of the spectrum is the existence of an S = 1 two-magnon bound state 
for a range of q around q ~ (2n + l)ir where n is an integer. This bound state is due to 
the attraction between adjacent excited dimers with total spin 1 (and 0). It exists only 
over a limited range of q around the bandwidth minimum in the two-magnon continuum 
because of the nature of the hopping matrix elements. Near the two-magnon bandwidth 
minimum the hopping element is small, and binding occurs. Far from this minimum the 
gain in hopping energy of excited dimers is dominant and no S = 1 bound mode exists. 
An interesting consequence of the competition is that an S = bound mode (not visible 
to neutron scattering, but can be observed using Raman scattering) should exist below 
the continuum for all q, because it has a larger attractive interaction between adjacent 
dimers. 

The strength of production of the S = 1 bound mode in neutron scattering depends 
on the spin operator matrix element to this state. Perturbation analysis in a || shows 
that there is no zeroth order coupling of the bare dimer ground state to the two-magnon 
bound state and continuum through the neutron scattering matrix element. The leading 
perturbative contribution to S(Q,u) appears at 0(a 2 ), and is due to a transition from 
the 0(a) two-excited-dimer component of the ground state and to the 0(a) one-excited- 
dimer component of the bound state. The matrix element of these two basis components 
has a complicated and characteristic wavevector dependence. 

Figure (2) shows the energy corresponding to a weighted average of scattering as a 
thick grey line. It is notable that at q = 3tt to a good approximation the neutrons 
couple only to the bound mode, so that nearly all the scattering weight is in it, not the 
continuum. The calculated neutron scattering intensity from the bound state is 2% of the 
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one-magnon intensity which agrees with the data in Figure (1). 

The one- and two-magnon scattering at 300 mK was scanned from q = it to 5tt in steps 
of 7r/4. The background subtracted data are plotted in the upper panel of Figure (3). 
A calculation of the magnetic scattering based on the perturbation calculation described 
above, including the dimer envelope function is shown in the lower panel of Figure 
(3). The calculation is directly comparable with the data in the upper panel of Figure 
(3), and at a qualitative level there is good agreement with experiment. 

A quantitative comparison between theory and data is shown in Figure (4). The 
measured positions of one- and two-magnon peaks are plotted in the left panel. Energies, 
widths and intensities for each peak were extracted by least-squares fitting of Gaussians. 
Considerable dispersion of the one-magnon modes is evident. Measurements of the one 
magnon dispersion were previously fitted using an 0{a) model of chains with an additional 
weak interdimer coupling JL / 200 i and J[q 01 / 2 ] and . 



The dispersion relation predicted by this model gives a good account of our data, Fig- 
ure (4). The two-magnon peak is broader than experimental resolution, and the extracted 
positions (grey filled circles) are nearly dispersionless and are located at the calculated 
weighted average energies (grey band) replotted from Figure (2). Interchain effects are 
effectively integrated over in the two-magnon scattering, and this results in line broad- 
ening rather than shifts in energy; we thus have implicitly included interchain coupling 
effects in our definition of a = 0.27 for (1). 

The right-hand panels of Figure (4) show one and two magnon integrated intensities 
extracted from the fit. The one-magnon intensity (lower panel) is well described by the 
theory (solid line), and shows a simple oscillation. Because the one-magnon intensity 
is modulated by the dimer envelope function in S(Q,u) the intensity should go to 
zero near q = 4n. The residual intensity comes from secondary elastic scattering from 
incoherent processes, and the theory (solid line) is corrected for this. The two-magnon 
intensity (upper panel) shows a more complicated q dependence, which is due to the 
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various basis transitions that contribute to the coupling of the ground to excited states. 
Although the 0(a) result looks qualitatively similar to the data, it underestimates the 
scattering at q ~ |7r and overestimates it at 2n, which may indicate that higher order 
terms in the scattering amplitude are important. It is notable that the two-magnon 
intensity is very strongly dependent on the spatial arrangement of magnetic ions. 

The binding energy of the S = 1 state is predicted to be || E B = J — |§a 2 ) = 
0.017 meV for CN. Scattering around q = 37r is centered at 0.852 ± .007 meV, which gives 
Eb = 0.03 ± .02 meV. Although the energy and intensity around q = 3ir lend support 
to binding around this bandwidth minimum, the experimental error means this does not 
constitute definitive proof of the effect in CN. The much better energy discrimination 
of time-of-flight (TOF) neutron spectrometers could provide this by resolving the bound 
mode from the continuum. Previously, TOF techniques have proven successful in the 
study of similar binding effects at the bandwidth minimum of the two-soliton continuum 



scattering of the S = 1/2 XXZ Ising chain material CsCoCLj |T2|]. Unlike the AHO 
binding only occurs when extra terms in the Hamiltonian, such as exchange mixing W2 
or next-nearest neighbor coupling [I3|, are present. However limited neutron fluxes may 
make such measurements difficult for CN ||. 

In conclusion, we have used inelastic neutron scattering to investigate the ground and 
excited states of the near-ideal alternating Heisenberg chain material Cu(N03)2 ■ 2.5D 2 0. 
Our measurements are consistent with predictions of this model for several magnetic 
properties of this system, including the ground state energy, one- and two-magnon exci- 
tation spectra and intensities, and possibly the existence of a two-magnon bound state. 
Much experimental work remains to be done to establish the phenomenology of binding 
in isotropic ID systems. 
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FIG. 1. (a) Low temperature scattering at q = 2ir. The dashed line is a fitted background 
and the solid line is a fit to the scattering described in the text, (b) Two magnon scattering 
with background subtracted off. The solid line is a fit (see text). The solid bar indicates the 
instrumental resolution, (c) Two-magnon scattering for q = 3ir with nonmagnetic background 
subtracted off. 

FIG. 2. Schematic showing predictions of perturbation theory specialized to J = 0.45 meV 
and a = 0.27. 

FIG. 3. (Color) Upper panel shows a color filled contour plot of the measured data with 
nonmagnetic background subtracted. Intensity is on a linear scale indicated by color, going 
from dark red (minimum) to light yellow (maximum). The two magnon scattering has been 
enhanced by a factor of 10 2 to make it visible on the same scale. Lower panel shows the 
calculated scattering using perturbation theory with corrections for instrumental resolution, 
multiple scattering and magnetic form factor. 

FIG. 4. Comparison of theory and data. Left panel shows fitted positions of observed 
scattering with predictions using perturbation theory (see text). Right lower panel shows fitted 
one magnon intensity compared with perturbation theory (see text). Right upper panel shows 
a comparison of two-magnon intensity with perturbation theory (see text). 
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